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ABSTRACT 

A model for the coherent output coupler of the Bose-Einstein condensed àtoms from a trap in 
the recent MIT experiment (Phys. Rev. Lett., 78(1997)582) is established with a simple many- 
boson system of two states with linear coupling. Its exact solution for the many-body problem 
shows a factorization of dynamical evolution process, i.e., the wave function initially prepared 
in a direct product of a vacuum state and a coherent state remains in a direct product of two 
coherent states at any instance in the evolution of the total system. This conclusion always holds 
even for a system with a finite average particle number in the initial state. Its thermodynamical 
limit can be directly dealt with in the Bogoliubov approximation and manifests that an ideal 
condensate in the trap will remain in a coherent state after the r.f. interaction while the output- 
coupler pulse of àtoms is also in a coherent state, which means a coherent output of atòmic beam 
to form a macroscopic quantum state in a propagating mode. 
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1. Introduction 

According to de Broglie matter wave theory, all matter in general can behave itself like a light 
that spreads out in space and combines with another one to form interference pattern. In 
principle, it is logical to image that a "matter wave làser" can be realized, in analog to a light 
làser, to produce the output coherent matter wave . Based on such fundamental consideration 
and stimulated by the experiments observing atòmic Boson-Einstein condensate (BEC) in last 
two years [1-3], a number of efforts coming from different research groups have been placed 
on both its experimental setup and theoretical possibilities [4-7]. Until the end of last year, a 
milestone experiment for the matter wave làser of bosonic àtoms (also called atom làser) was 
accomplished by Mewes et al in MIT [8]. 

In their experiment, they coupled the toms trapped in a state |1) with that untrapped in 
another state |2) by a sweeping radiation frequency (r.f.) field. The r.f. field coherently turns 
the BEC àtoms constrained in |1) into a propagating mode of àtoms in |2). The experiment 
shows that " the condensate which is initially in a coherent state remains in such a state after 
the r.f. interaction while the output-coupler pulse of atom is also in a coherent state" [8]. In 
fact, a làser of photons above threshold can be understood in terms of the coherent state for its 
well-defined classical analogy with minimum quantum fluctuation in amplitude [9]. Similarly, it 
is reasonable to believe that such a output of àtoms in a single coherent state is something like 
the photon làser. Roughly, we can regard the dynamic process of wave function evolution driven 
by the r.f. field as a "stimulated" process of (massive) matter wave. 

In this paper, we first use the Bogoliubov approximation [10,11] in thermodynamic limit to 
give a theoretical analysis for the dynamic lasering process of BEC àtoms in the MIT experi- 
ment. We will show that such an approximation is not actually needed to reach the conclusion 
about coherent output in the requirement of very large N c . Namely, even for the case that the 
number N c of àtoms condensed in the ground state is not very large, thereby the Bogoliubov 
approximation is violated, the output of finite àtoms is still coherent so long as the system is 
interatomic interaction free and the àtoms in the trap is initially in a coherent state. Theoreti- 
cally, this is the factorization structure for the evolution of the many-boson model of two state 
with linear coupling, i.e., the wave function at any instance remains in a direct product of two 
coherent states if the system is initially prepared in a direct product of a vacuum state and a 
coherent state. This fact does not depend on the average number of condensed àtoms in the 
initial state and it implies that the experiment result by Mewes et al is quite steady even for 
rather finite N c . 

2. A Dynamic Model of Output Coupler via Bogoliubov Approximation 

The second quantization of model Hamiltonian for the MIT experiment can be written as 



in terms of the creation and annihilation operators, b{, b 2 , &i an d &2, of bosonic àtoms for the 
magnetically trapped state |1) and the untrapped state |2) with level difference hu a . The r.f. 
pulse of varying frequency u(t) is a sweeping classical electromagnètic field, which couples |1) 




(1) 



and |2) through the dipole matrix element ku R = Jhu/2e V = hg/\^V. V is the effective mode 
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volume and Eq is the vacuum permittivity. This many atom system without interatomic interac- 
tion can be modeled as a linear coupling system of two-oscillators. The main simplification is to 
ignore the quantized motion of atòmic center of mass in the trapped state by an inhomogeneous 
magnètic field. As pointed out by Mewes et al, this inhomogeneous feature can be neglected for 
a sufficiently short r.f. pulse. 

The MIT experiment is theoretically described as an initial value problem for Schòrdinger 
equation governed by H with the initial state |t/>(0)) = \a = V^c)i ® 1 0) 2 - Here, \a = V^Qi is a 
Glauber coherent state of the operator b\ characterizing N c àtoms condensed in the trapped state 
|1). No àtoms occupy the untrapped state |2) initially. Why does the coherent state \a) represent 
the BEC is still an open question, but there are some reasons enable one to believe it is right, 
such as the phase locking implied by BEC [12] and the correct average atòmic number N c . The 
additional reason we believe is due to the connection with the Bogoliubov approximation [10,11], 



which replaces operators b\ and òj with a c-number yN c when much many àtoms condensate in 
the trapped state. A factorizable evolution structure in the MIT experiment is just rooted in this 
connection and will be discussed as following. There is a close relation between the selection of 
initial coherent state \a) and the Bogoliubov approximation. In fact, a unitary transformation 

\(f)(t)) = D~ 1 (a)\ip(t)), in terms of the generator D(a) = expfaòj — abi] of coherent state, 
defines an equivalent initial value problem of Schòrdinger equation for |0(O)) = |0)i ® |0) 2 and 
the equivalent Hamiltonian 



D(a)- 1 HD(a) =H + hto R J N ( 



b\ exp ^— i J üj{r)dT^j + h. 



For the case of very large N c in BEC, the coupling term is very small in comparison with the 
term of fiuiR\fW c and can be neglected to obtain an effective Hamiltonian 



Hb = hu a b\b 2 + hujR\lN t 



b\ exp ( — i I üj{r)dT 



+ h.c 



(2) 



This is just the Bogoliubov approximation of the original Hamiltonian H. 

Now let's consider the thermodynamical limit of an infinite number of àtoms in an infmite 
volume but with the density fixed, i.e., n c = \im N y^oo (N c /V) — > constant, the coupling term 
hujR oc 1/VV — > 0, but hujR\fW~ c = Tig^/n^is finite. Therefore, in the point of view of Schòrdinger 
evolution, the Bogoliubov approximation for our problem just describes an excitation motion of 
system in a large background provided by the BEC initial state. Namely, for a very large iV c , 
the system initially in a coherent state remains in such a state while another component \&{t)) 
is governed by the Bogoliubov approximate Hamiltonian. Mathematically, the evolution wave 
function enjoys a factorizable structure 



m)) = w)®m))- 



(3) 



The similar factorization structure has been shown in the analysis about the transition from 
quantum domain to macroscopic, classical domain concerning the wave function collapse in 
quantum measurement [13] and the macroscopic damping in quantum dissipation [14,15]. The 
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method calculating wave function developed in [14,15] will be used in this paper to obtain the 
exact solution for factorized wave function in section 4. 
3. Coherent Output via Factorization 

The next step is to prové | <£>(£)) is still a coherent state evolving from its vacuum state 1 0)2- 
It is quite clear that the Bogoliubov approximation Hamiltonian (2) is just a forced harmònic 
oscillator (FHO). It is well-known that the displacement property of a quantized FHO results 
dynamically in a coherent state for the initial vacuum state. By considering the case of uj(t) 
independent of time, i.e., uj(t) = ou, the solution of the Hamiltonian (2) can be found in some 
Standard textbook [16] = \à(t)) with 

& (t) = ^5 ( e >A t _ ^ e-toj (4) 

Namely, the total wave function in the BEC with large iV c can be factorized in the direct product 
of two coherent states \ip(t)) = \a) ® |<5(í)). The second components implies the output-couple 
pulse of the atòmic beam is in the coherent state, which "can be regarded as a pulsed atòmic 
làser" . 

Notice again that the above theoretical discussion about the output- coupler for atòmic làser 
is only with the help of a highly-simplified many- boson model of two states with linear coupling 
rather than the complicated nonlinear one, such as the system of Gross-Petaevskii equation. 
Because the atòmic number of the coherent output in the untrapped state 

(mibWim) = \à\ 2 = 2^fiVc(l -cosAí) (5) 

is an oscillating function of t, the sweeping mechanism [17, 8] of the r.f. pulse must be invoked 
here to localized the population in the untrapped state at t = 00. This will be in the aid 
of Landau-Zener method [17] to control the populations in trapped and untrapped states by 
adjusting the time- dependent frequency u(t) of the r.f. pulse from diabètic to adiabatic point 
in the avoid level crossing. A direct treatment for the r.f. sweeping is changing the factors 
exp (±iut) into exp (±.i J* * üj(t)<1t^ of H in the above approach. Near the resonance uj(t ) = cü a , 

A~o;(í )(í-ío), 




This leads to a fixed coherent population in the untrapped state, which can be controlled by 
adjustment of the changing rate òü(t ) of the frequency of the r.f. pulse at resonance. 
4. Coherence From Finite Condensed Àtoms 

The above discussion shows that the factorized structure of wave function crucially result in the 
output of coherent atòmic beam as an atòmic làser when the initial state of system is prepared 
on the BEC in a trap for which the Bogoliubov approximation is hold. This fact seems that 
the factorizable structure of evolution depends on the Bogoliubov approximation, namely, it is 
decided by whether or not there are very many àtoms in the BEC trapped state. In the following, 
it will be demonstrated by the exact solution of the Hamiltonian (1) that the requirement of 
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very large N c is not necessary so long as the àtoms are prepared in a coherent state \a) even 

with rather finite average atom number N c = {a\b\b2\ct) . The coherent state is only a special 
Poisson superposition of infinite number states with marched phases 6 n = n6, n — 1,2, ... . 

The second reason to consider the exact solution concerns the analogy in the photon case. In 
the theory of light làser, a sufficient strong light field can be treated as a classical electromagnètic 
field and it can be regarded as a coherent state above the threshold [9]. In this sense, the quantum 
effect mainly results from the quantum fluctuation, which is demonstrated by a weak light field. 
Therefore, for the matter wave of massive partides, the infinite number case only corresponds 
to the classical analogy of photon field. Then, we need to study the case corresponding to the 
quantum fluctuation of the weak light field, for which the Bogoliubov approximation is broken 
down. 

According to [14,15], the exact solution of the wave function of the linear system in the 
Schròdinger picture can be easily obtained from that of the canonical operators, such as 6i(í) 
and 6 2 (i), in Heisenberg picture. In the case of ui{r) independent of time, by solving a system of 
one-order differential equations resulting from the Heisenberg equations of bi(t) and b 2 (t), the 
exact expressions of these solutions are 



where 



Pi(t) 

M) = 

ai(í) 



6i(í) =a 1 (t)b 1 (0) + a 2 (t)b 2 (0), 
h(t) =/3i(í)6i(0)+/3 2 (í)6 2 (0) 



a 2 (t) = 



^A 2 + Au 

1 r 
A 2 + 44 
1 

\M 2 + 44 



■iuit 



-ÍUlt 
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^A 2 + Auj 



(7) 
(8) 

(9) 
(10) 

(11) 
(12) 

(13) 



From the property of the evolution operator, U (t), of the total system,i.e., b a (t) = U (t)^b a (0)U(t), 
(a = 1,2), the initial state \ip(0)) — \ct — \fW c ) <8> |0) will evolve into 



i^(í)> = u(t)\m) = e-^ 2 e zjbii-m ® io> 

n=0 a ' 

00 n / t 4. \ n 

= E - ( ai (-m(0) + a 2 (-t)bl(0)) |0) ® |0). 



(14) 



n=0 
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That is just a direct product of two coherent states 

\ip{t)) = |ac*i(-í)) <g> \aa 2 (-t)) 

= | jN~ c e-ï At (cos(ííí) + i sin(fií) cos 0)) (g) | - iy^e - ^ Aí sin sm(Qí)e" íwí ) (15) 

where = 2cü r /A, where fi = y / A 2 /4 + o;|. 

This result indeed leads to a coherent output of àtoms in the untrapped state so long as the 
àtoms are initially prepared as the BEC in trapped state even for the finite average number of 
àtoms. In this output state with wave vector k, the expectation value 

(0( x )) = (^(í)|0(x)|^(í)) = J ^ sin 9 sin(í)í) sin (kx - ^ At - wí) (16) 



of the bosonic field, 0(rr) = J2k \Jh/2uV (b 2 ke lkx + h.cj, looks just like a classical coherent field, 

where the variance of the field A0 = \J (0(x) 2 ) — (4>(x)) do not depend on N c . This means the 
relative uncertainty of the amplitude, A0/ (<f)(x)), will approach zero only in the thermodynamics 
limit: N c , V — > oo, but N c /V — > constant. Therefore, for finite iV c , the quantum fluctuation 
A(f) oc \fh can be neglected for the dynamic problem. The Bogoliubov approximation in the 
last section is just the thermodynamics limit of the results in this section. Notice that, in the 
thermodynamics limit, we have ojr — > 0, Vt — > y, but N c sm.9 oc do not approach zero. 

5. Discussions 

In summary, for the MIT experiment on the output coupler of atòmic làser, we have firstly shown 
that, if the àtoms with an average number N c are accumulated in a single quantum state |1) to 
form a coherent state \a = ^/N~ c )i and no àtoms occupy the untrapped state |2), the spin-like 
dynamics of Rabi oscillation will create a factorized wave function \ip(t)) = |a(í))i<g) \(3(t)) 2 , which 
is a product of two coherent states. Even for the case without ideal Bose-Einstein condensation 
that N c is not sufliciently large and thus \a = V^c)i is n °t a macroscopic quantum state, 
this wave-function of evolution remains in a factorizable structure so long as \a = V^c)i is 
coherent states. However, it seems very difficult to prepare finite numbers of àtoms in a coherent 
state duo to the massive feature of àtoms. The spin-like dynamics of Rabi oscillation will only 
create an entanglement evolution state | ?/>(£)) = \ a j(t))i ® \Pi(t))2 for an arbitrary initial 
components |^(0))i other than a coherent state. Only in BEC case that the fixed phase and 
amplitude of macroscopic wave function are defined and thus to form a Glauber coherent state, 
this entanglement state can be completely factorized. 

Such a factorized structure of wave function first pointed out by Mewes et al is very crucial 
to realize their experiment, but it dynamical origins mainly depend on the linearity of coupling. 
For the case with an interatomic interaction in the trapped or untrapped states, which roughly 

is of the form bjbjbibi, the corresponding Heisenberg equations are no longer linear and thus this 
factorized structure will be broken down. Therefore, the sufliciently weak interaction should be 
required for a coherent output in practical system. In fact, Ballagh et al [18] have modeled a 
theory for the evolution of two component BEC by generalizing the Gross-Pitaevskii equation. 
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Their theory introduces a non-linear element resulting from the interatomic interaction and 
shows a rich variety of phenomena of nonlinear effect such as the localization of atom in the 
untrapped state due to the nonlinear coherent coupling. 

Notice that an interference experiment was also carried out by Andrews et al in the same 
research group at MIT [19] to examine the coherence of àtoms in BEC by overlapping two 
expanding condensates flying from the double well magnetically-trapped potential when this 
trap is adiabatically switched off. The next step following this experiment was immediately 
accomplished to test the coherence of the output couplers rather than a pair of expanding atòmic 
sources. They also have shown that a similar high-contrast regular pattern of interference. Based 
on an intuitional understanding of the MIT "atòmic làser" in terms of coherent states described 
in this paper, a many-body quantum theory for the interference of output couplers can be 
present in a forecoming paper as a space-dependent generalization of present results. The non- 
vanishing order parameter of quantum system, concerning a symmetry broken involving the 
atòmic number conservation, can be obtained for the description of regular pattern of stripes in 
these interference experiments. 

One (CPS) of t he authors wish to express his sincere thanks to K. Young for inviting him 
to visit the Chinese University of Hong Kong as a C. N. Yang's Fellow. The work is supported 
in part by the NSF of China. 
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